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We investigate the energy of a theory with a unit vector field (the "aether") coupled to gravity. 
Both the Weinberg and Einstein type energy-momentum pseudotensors are employed. In the lin- 
earized theory we find expressions for the energy density of the 5 wave modes. The requirement that 
the modes have positive energy is then used to constrain the theory. In the fully non-linear theory 
we compute the total energy of an asymptotically flat spacetime. The resulting energy expression 
is modified by the presence of the aether due to the non-zero value of the unit vector at inflnity and 
its 1/r falloff. The question of non-linear energy positivity is also discussed, but not resolved. 



a^ ■ 
o : 
o , 

(N . 

o ; 

o- 

> 

o\ : 

in . 
o 
i> 
o 

vn . 
o . 

"o ■ 
cr 

u 



I. INTRODUCTION 

Although Lorentz invariance has been a key feature of 
theoretical physics for a century, recently there have been 
a number of reasons for questioning whether it holds at 
all energy scales. For example, some possible quantum 
gravity effects hint that it may not be a fundamental 
symmetry Thus, it is useful to construct effective, 
low energy symmetry breaking models in the regimes of 
the Standard Model and General Relativity (GR). The 
new effects that appear can then be studied in a fa- 
miliar context and compared with observations. In the 
fiat spacetime background used in the Standard Model, 
Lorentz invariance can be broken by background tensor 
fields [2] . However, when we attempt to couple such fields 
to gravity, they will also break the general covariance of 
GR, which we regard as fundamental. In order to by- 
pass this problem, it is straightforward to consider these 
fields as dynamical quantities along with the metric. In 
this paper, the source of the Lorentz violation (LV) will 
be modeled as a unit timelike vector u". The unit time- 
like restriction preserves the well-tested 5*0(3) group of 
rotations while enforcing the breaking of boost symmetry 
at every point in the curved spacetime. Therefore, the 
vector can be said to act as an "aether" . 

Similar work was initiated in the early 1970's by Will, 
Nordtvedt and Hellings [1, 0, \3\ who studied a vector- 
tensor model without the constraint in the context of al- 
ternative theories of gravity. For a review of more recent 
work on the subject, including aspects of observational 
constraints, waves, cosmology, and black holes, see @ 
and the references therein. Following these authors, we 
will refer to this theory as the "Einstein- Aether" theory. 
One important open question is whether the Einstein- 
Aether theory is energetically viable and stable. The 
Will-Nordtvedt-Hellings models, for example, are unsta- 
ble because fluctuations of the unconstrained vector can 
be either timelike or spacelike, allowing ghost configura- 
tions and energy of arbitrary sign [2] . 

Energy in a field theory is defined as the the value 
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of the Hamiltonian, which acts as the generator of time 
translations. Although in diffeomorphism invariant the- 
ories there is generally no preferred notion of time (and 
thus energy), in asymptotically flat spacetimes one can 
naturally define the ADM and Bondi energies associated 
with asymptotic time translations at spatial and null in- 
finity respectively. The ADM and Bondi definitions for 
GR have also been shown to satisfy positive energy the- 
orems [1]. 

In this paper we examine energy in the Einstein- Aether 
theory. Since it has proven difficult to directly construct 
the Hamiltonian for the theory, we instead consider the 
pseudotensor method of studying gravitational energy. 
Such an approach was first taken in a similiar context by 
Lee, Lightman, and Ni who derived pseudotensors for 
the unconstrained vector-tensor models but did not eval- 
uate them on solutions. Despite the non-covariance of 
pseudotensors, it is known that they give well-defined re- 
sults for the spatially averaged energy carried by waves in 
linearized theory and the total energy of asymptotically 
fiat spacetimes. In gravitational wave physics they pro- 
vide a simple and straightforward method for calculating 
averaged energy densities and the energy-momentum flux 
radiated away from sources. In addition, Ghang, Nester, 
and Chen [l^ have shown that the superpotential asso- 
ciated with every pseudotensor corresponds to a (albeit 
non-covariant) quasi-local Hamiltonian boundary term. 

We first discuss the Einstein- Aether theory and then 
motivate and construct its modified Weinberg pseudoten- 
sor expression. As the calculational and consistency 
check we also use a Lagrangian based method to de- 
rive the modified Einstein "canonical" pseudotensor and 
its associated superpotential. We then apply these ex- 
pressions to solutions in both the linear and non-linear 
regimes. In the linearized theory we find that the Ein- 
stein and Weinberg prescriptions give the same energy 
densities for the plane wave modes derived in [ifj. Re- 
stricting these densities to be positive yields constraints 
on the model in terms of the coefficients of the aether part 
of the action. These constraints are also compared to re- 
sults obtained in the limit where the metric and aether 
decouple [l2|- In the full non- linear theory the Einstein- 
Aether superpotential is used to obtain the total energy 
for an asymptotically fiat spacetime. This result agrees 
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with fl§\, where the total energy in the Einstein- Aether 
theory is derived via the covariant Noether charge for- 
mahsm. We conclude with a discussion of the status of 
positive energy in the non-linear regime and prospects 
for a positive energy theorem. 



II. EINSTEIN-AETHER ACTION 



^ ab = Vm[J{a Uf,) - J (^a^b) - J (ab)U ) 

+Cl [(V„U„)(V'"ub) - (V„7i,„)(Vbli™)] 
-I-C4 ilaUb 

+ [u«(V™ J""') - ail^] UaUb 

1 ^ ^ 

-7;Lu9ab, (9) 



We can model gravity with a dynamical preferred 
frame using a timelike unit vector u". This vector field 
breaks local Lorentz invariance spontaneously in every 
configuration leaving behind the 3-D rotation group as 
the residual symmetry. The unit norm condition is re- 
quired to avoid ghosts and also desirable because we re- 
gard the norm as extra information beyond what is nec- 
essary to determine the preferred frame. Taking an ef- 
fective field theory point of view we can consider the ac- 
tion as a derivative expansion, subject to diffeomorphism 
symmetry. The result, up to two or fewer derivatives, is 
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where 



-R - if^^^^^VaW^Vfeu" - Xigabu'^u' - 1). (2) 



The "kinetic" term K"-''mn is defined as 

if^'mn = cig'^'^gmn + CsC'^' + CsS^^S^ + du'^u^'g^^ (3) 

and A is a Lagrange multiplier enforcing the unit timelike 
constraint. R is the familiar Ricci scalar and the coef- 
ficients Cj in K'^^mn are dimensionless constants. Note 
that a term proportional to RabU°'u^ is not explicitly in- 
cluded as it comes about as a combination of the C2 and 

C3 terms in ([1]). The metric signature is (-1 ) and 

the units are chosen so that the speed of light defined by 
the metric gab is unity. 

The field equations from varying the action in ([T]) to- 
gether with a matter action with respect to g"^ and u° 
are given by 

Gab = T^at^ + SnGT^'ab (4) 

^ aJ'^m ~ C4UaVmU'^ — Ali™, (5) 

(6) 



gabu'^u'' = 1. 



where 



and 



<^ m — ^ ran v ci 



Ua=U VbUa- 



(7) 
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Here we assume that there are no aether-matter couplings 
in the matter action. The aether stress tensor is given by 



Lagrange multiplier A has been eliminated from ^ by 
solving for it via the contraction of the aether field ([5]) 
with u°-. As we will see below, the the form of the aether 
stress tensor and Einstein- Aether Lagrangian will be im- 
portant tools in derivation of the modified Weinberg and 
Einstein pseudotensors. 



III. WEINBERG PSEUDOTENSOR 

Weinberg's pseudotensor construction is based on 
the "field theoretic" approach to GR that treats gravity 
as a spin-2 field on a flat background spacetime. Using 
Greek indices to represent coordinate indices, we begin by 
writing the metric in coordinates such that g^i, = rj^^i, + 
/i^jy, where 77^^ is the flat Minkowski metric and /i^^ is 
an symmetric tensor field with the asymptotic conditions 
V ~ 0(l/r), d^h^, ~ 0(l/r2), ~ 0{1/t^). 

The Einstein tensor can be expanded into a series of parts 
linear, quadratic, and higher order in the field variable 
/w. Following Ch. 20 of Misner, Thorne, and Wheeler 
[Iq the non-linear corrections to the Einstein tensor are 
defined as follows 



16^G V = 2G«-2G^„ 



(10) 



where G^^} and G^^ are the linearized and full non-linear 
Einstein tensors respectively. Note that this splitting is 
non-unique because it depends on the coordinate system. 
Since the linearized Einstein tensor is symmetric and sat- 
isfies a linearized Bianchi identity 9^G^V^ = 0, it can be 
rewritten in superpotential form 



2G(^) 



,a/3 



(11) 



where H^j^^ap has the symmetries of the Riemann ten- 
sor Hf^^aii = Hypii,][af3] = Hap^v (see, for example (IJl). 
Using pO|l and (fTTj) the full Einstein equation becomes 

i/^a./s'"'^ = 167rG {t^, + T^,). (12) 

Due to the symmetries of H^avf}) this implies that 
d^'it^.y + T^^) = 0. Therefore the integral of too + Too 
over a spacelike slice is a conserved quantity. This con- 
served quantity 



(too + Too) d^x 



1 



167rG 
1 

IGttG 



i?oao/3'"n'' d?x, (13) 
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where is the unit normal to the surface at spatial in- 
finity, is in fact the total energy, with tpo acting as the en- 
ergy density of the gravitational field alone. To sharpen 
this point, consider the case where the gravitational field 
is weak everywhere, allowing use of the linearized theory. 
The leftmost member of ([T3|) then gives the total mat- 
ter energy, which in this case is the total energy. The 
rightmost member is insensitive to the interior volume, 
so replacement by arbitrary sources and strong fields in 
the interior will not affect the identification of as the 
total energy. 

The extension to the Einstein-Aether theory is 
straightforward. The metric field equations (|4|) take the 
form 



G^. = G^, - ni;' - SttGT^,. (14) 



In addition to the metric, we now decompose the aether 
into background and dynamical part by writing u'' — 
_|_ Unlike normal matter fields the aether stress 

(u) 

Tfj^u contains linear pieces in the perturbation due 
to the fact that the aether does not vanish in the back- 
ground (since it is always a unit vector). These linear 
terms will modify the Weinberg pseudotensor and super- 
potential. Performing the split of the modified Einstein 
tensor Gab as in pU]) we find 



16^G V = "^G^^u 



2Gui/ 



(15) 



where G^i^ = G^^^-^ - ^(iX")^^. G^^ satisfies a Bianchi 
identity S/^Gf^i, if the aether is uncoupled to the 
matter and if the aether field equation ([5]) is satisfied. 
Therefore in the linearized case we can write 



2&' 



H 



,a(3 



along with 



,al3 _ 



167rG(V + 7^M.)- 



(16) 



(17) 



By the same reasoning as before we could conclude that 
the total energy is given by 



E 



1 



16nG 



(18) 



However, unlike the GR case (IT^ . it is not clear whether 
the new Weinberg superpotential H^aup'" can be ex- 
pressed as a local function of the fields hab and u° [l^l ■ 
On the other hand, the pseudotensor can be calcu- 

(u) 

lated directly via the non-linear pieces of T/iiy and G^^ 

in Gfj,u- This will be used to compute the linearized wave 
energy densities. Evaluation of the total energy as a sur- 
face integral at spatial infinity requires a locally defined 
superpotential. Since we do not have knowledge of the 
aether corrections to Weinberg superpotential we shall 
instead consider the Einstein superpotential, which can 
be derived directly from the form of the Lagrangian. The 
Einstein formulation of gravitational energy-momentum 
will also provide a consistency check when we evaluate 
the energy density of the linearized plane wave modes. 



IV. EINSTEIN "CANONICAL" 
PSEUDOTENSOR 

The gravitational energy pseudotensor originally de- 
rived by Einstein in 1916 shortly after his discovery of 
the field equations of GR is closely related to the familiar 
canonical stress tensor of matter fields in flat spacetime. 
In order to derive the corresponding expression for the 
Einstein- Aether theory, we use a Lagrangian approach 
based upon the famous work of Noether relating symme- 
tries to conservation laws. In flat spacetime, invariance 
of a Lagrangian under global space and time translations 
is associated with the conservation of energy-momentum 
expressed by the conservation of the canonical stress ten- 
sor 



rp fl 



dC 



(19) 



where C — C{il}, d^j) and ip represents a general collec- 
tion of flelds with indices suppressed. In the case of local 
symmetries, such as the diffeomorphism invariance of the 
Einstein- Aether theory, the situation is more complex. In 
the Appendix we review a general formalism due to Julia 
and Silva [l^ for constructing Noether currents and su- 
perpotentials and apply it to the Einstein- Aether theory. 
Using these results we show that the pseudotensor and 
superpotential have the following general form 



f ^ 
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dL 



IQttG \d{df,gai3) 
dL 



f—g( dL 



167rG \d{df^gai3) 
dL 

-o„.u ' 



5^L 



(20) 



(21) 



where L is the Lagrangian 

-K'^^^^Vo^u'^Vpu" - X{g^,,u^'u'' - 1). (22) 



T — — «"'Vr'' T" — r'' r"^ 



Note that we have eliminated a surface term in the 
Einstein-Hilbert Lagrangian, replacing the Ricci scalar 
R with the Einstein-Schrodinger 'T^" action, which de- 
pends only on the metric and its first derivatives. When 
evaluated on-shell the pseudotensor and superpotential 
obey the following relations 



(23) 
(24) 



To account for the presence of any non-aether matter 
sources one only has to make the replacement t^^ 

+ T^,^' in and ([Ml). Like the Weinberg construc- 
tion, the pseudotensor t^^ is a conserved quantity and is 
related to the divergence of a superpotential. 



4 



The contributions from the pure GR Lagrangian 

are the Einstein pseudotensor 



_i_r/3 p7 „Ma _ -pP y-y „Ma i -pa p/J f3y 



I T^p p/3 _ r,pp pa /37 

pa^ P7y 



and the von Freud superpotential (see, e.g. |19| ) 



(25) 



irrr A. 



^JPrd,{{-g){g^^g^-'-g^^g^'')}. 



(26) 

To compute the additional aether modifications, we use 
the relation 

d{d,g^0)-'^^ ^-^^'^^ ^'^-^ ^ ^'^9(r;,) 

(27) 

and L„ = K"''^ au'^V fju" in (HOI) and ^ since when 
evaluated on solutions any terms related to the unit con- 
straint will vanish. We find the pseudotensor 



1 



^gJ^pV^uP 



-g{(J^p + Jp^)u- 



-{J" (3 + J^u^ + (J^" - J"^)w/3}r\. 

+5^y^L„) . (28) 



and the superpotential 
1 



-{J'^p + J/3")w^ + (J^" - .r^)up) (29) 

where J°/3 is defined in ([7]). The above decompositions 
of (|20p and (|2ip into GR and aether pieces do not satisfy 
([23|) and ([24|) independently. A key requirement when 
evaluating these pseudotensorial expressions is that the 
metric must be written in a coordinate system where the 
connection coefiicients vanish like 0(l/r) or faster in the 
asymptotic limit. If the coordinate system is not chosen 
properly then these expressions will yield incorrect ener- 
gies and momenta j25j . This condition was not well un- 
derstood in the early literature on gravitational energy- 
momentum, but can now be explained using an analysis 
of the boundary terms and conditions in an action. See 
the Appendix for further details. 



V. ENERGY IN LINEARIZED THEORY 

Equipped with the modified Einstein and Weinberg 
pseudotensors we can now calculate the energy density 
of the linearized plane wave solutions to the Einstein- 
Aether theory. The plane wave solutions in the absence of 
matter are found by linearizing the field equations above. 



gl)-®, with .g^^ 
gives 



hfj_i, and 



(1) 



,,0 



Q/3 

1, 



This 

(30) 
(31) 
(32) 



Cartesian coordinates are used in the flat background, 
r]^^ = (1,-1,-1,-1) and w'' = (1,0,0,0). Since the 
background value of the Lagrange multiplier vanishes, A 
in (|30p represents a perturbation. The superscript (1) 
represents quantities written to first order in the per- 
turbation. Jacobson and Mattingly [ll| then proceed to 
analyze these equations using the gauge choice 



ho., = 

ViA = 



(33) 
(34) 



which they prove to be accessible. Inserting plane wave 
solutions 



(35) 
(36) 



into the equations of motion, imposing the 4 gauge con- 
ditions p3|) - (p4)) . and choosing coordinates such that the 
wave- vector is (fco, 0, 0, k^) (travelling in the z direction), 
it is found [ll] that the mode polarizations and speeds 
are completely determined. The result is a total of 5 
wave modes falling into spin 2, spin 1, and spin types 
as shown in Table |T] [2^ . The notation / in subscript 
refers to the transverse components of the metric and 
aether while C14 — ci + Ci, etc. The 2 spin-2 TT metric 
modes look exactly like the usual GR case, except for the 
modification of the speed. The 2 spin-1 transverse aether 
modes and 1 spin-0 trace mode are new modes coming 
from the constrained aether, which is characterized by 3 
degrees of freedom. 

In order to determine the energy, note that in the ab- 
sence of matter the Weinberg prescription reduces 
to 



E : 



1 



16nG 



HoaOP 



inn cPx, 



(37) 



which clearly produces infinite total energy for plane 
wave modes. One could reformulate the problem in terms 
of wavepackets with the appropriate asymptotic fall-off 
conditions, but a far more direct approach is to simply 
evaluate the plane wave energy density too. This quantity 
is meaningless at a point for plane waves, but the aver- 
age over a cycle is well-defined. Consider a large, but 
finite region with nearly plane waves. There are "surface 
effects", but the the contribution to J tQ^d^x is domi- 
nated by the volume. Thus, too gives an effective energy 
density. 

The 3 general classes of modes were analyzed sepa- 
rately using the Riemann tensor package [2d| in Maple. 
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TABLE I: Wave Mode Speeds and Polarizations 

Mode Squared Speed Polarizations 

spin-2 1/(1 -C13) hi2,hii^h22 

Spin-1 (ci - ic? + ic3)/ci4(l - C13) hl3 = [Ci3/{1 - Cl3)s]vi 

spin-0 £123(2 - Ci4)/ci4(l - ci3)(2 + Ci3 + 3c2) hoo = -2vo, 

fill = ^22 — —C14V0, 

h:y.i = [2C14(C2 + 1)/Ci23]«0 



The package allows the user to enter the components of 
the metric and aether vector, calculate curvature tensors, 
and to define new tensors involving both ordinary and co- 
variant derivatives. In this case we entered the linearized 
metric and aether, where hf^^ and take the plane wave 
forms. A polarization was written as 

Aexpi/c3(z — st) + 74exp(— ifc3(z — st)) (38) 

where s are speeds shown in Table H] and A is a complex- 
valued function. Using this metric we calculated the ex- 
plicit form of the Weinberg pseudotensor ioo (flS]) up to 
quadratic order. Higher order terms will be small in the 

linearized theory and oscillatory terms proportional to 
2 

A and can be neglected in the usual time averaging 
process. These energy densities were then compared with 
the modified Einstein pseudotensor '''"''*'^'"tQ°+'^io° from 



(^5)1 and (^5)1 again up to quadratic order in the perturba- 
tions. Note that while ((23)l holds at quadratic order when 
the linearized equations of motion are imposed, (j24p does 
not. Therefore, one must use the modified Einstein pseu- 
dotensor directly to compute the energy densities. The 
results of the Weinberg and Einstein prescriptions agreed 
and are displayed below: 



^spin— 2 
^spin— 1 
^spin— 



SttG ^ ' ' 



1 



SttG ' ^ 
1 



kl \A\ 



2 "-s 



+ 2ci 



1 - Cl - C3 



SttG 



kl \A\^cu{2-cii) 



(39) 
(40) 
(41) 



These results have been independently verified in [2l| 
using the Noether charge method and a decomposition of 
u° into irreducible pieces. The lack of q dependence in 
([39|) and the simplicity of ([40 |l - (|4T|) is striking considering 
the complicated form of the pseudotensor expressions. 
The energy of the spin-2 mode is positive definite, like 
pure GR, while for the other 2 modes the sign of the 
energy density depends upon a combination of ci, C3, and 
C4. Note that when the q's are zero (HUl) and (|^T|l are 
zero as expected. This set of results for the coefficients 
also holds for exponentially growing modes (i.e. when 



<0) 



A cos(fcz + ip) exp{kst) 



(42) 



when we average over the spatial oscillations. Restricting 
> in Table |T] to eliminate the unstable modes and 



enforcing positivity in the energy densities in (I39p - (l4ip 
restricts the Ci values in the Einstein- Aether theory. 

In p^ . Lim worked in the limit where the aether and 
metric perturbations decouple, with the aether propa- 
gating in flat spacetime. Mathematically this amounts to 
tuning Ci, G — > while holding the ratio q/G fixed in the 
action ([1]). If we then expand the metric as g = 77 + VGh 
and take the limit, the action reduces to that of linearized 
gravity plus aether terms coupled only to r]ab- In this 



limit the linearized constraint reduces to 



= and 



we can decompose v'' into spin-0 and spin-1 parts via 
= d'^S + TV* where N\ = 0. By examining the Hamil- 
tonian of these modes they found ci > for positivity 
in both cases, neglecting C4. We can make contact with 
this result simply by examining in the small c; limit of 
the wave solutions. The trace and transverse aether en- 
ergy waves then correspond to the flat spacetime spin-0 
and spin-1 modes. To lowest order in Ci/G we find that 



Ci4 > 
Cl > 



(43) 
(44) 



for positive energy densities of the spin-0 and spin-1 
modes respectively. Restoring C4 in the flat spacetime 
analysis yields complete agreement. Note that for small 
Ci the > criteria for stable, non-exponentially grow- 
ing modes reduce to Ci/ci4 > for the spin 1 aether- 
metric mode and C123/C14 > for the spin trace mode. 
Thus, modes with positive energy are stable if C123 > 0. 



VI. NON-LINEAR ENERGY 

In this section we will attempt to extend the criteria 
for positive energy from linearized theory into the non- 
linear regime. As a first step, let us consider the total 
energy of an asymptotically flat spacetime in the full non- 
linear theory. Integrating ((24)) over a spacelike slice in the 
presence of non-aether matter gives the total energy 



cffO 



dx '°'Uo 



AO 



C/o "^^Ad5• (45) 



where '°*[/ = ^' U are the aether and von- Freud su- 

perpotentials, (^5)1 and (pS)) . and Tgg = t+T is total mat- 
ter and gravitational energy-momentum. The problem 
now is to calculate the superpotentials for the asymptoti- 
cally flat solutions to the Einstein- Aether theory. We will 
use Cartesian coordinates throughout since these have 
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the required asymptotic behavior discussed at the end 
of Section IIVI Therefore, the surface element is dS — 
r^dfi^ and the unit normal is (1/2, a^/r, y/r, z/r) where 
r = x-^ + + z'^ . For asymptotically flat boundary 
conditions we will assume that as r 



00 



Qtiy = fitly + 0{l/r) + 
u'^ = u^ + 0(l/r) + 



(46) 
(47) 



where y}^ = (1,0,0,0) with respect to the Minkowksi 
metric ry^,/ = (1,-1,— 1,-1). Equation (j45|) will only be 
affected by terms in the metric and aethe r up to 0(l/r). 
Using the analysis of the Newtonian limit [22l | and apply- 
ing the unit constraint, we find that far from the source 
in any asymptotically flat solution 



.900 

9i3 



-1 



r 



0(1/H 

2r 



1 



= 0(l/r2) + 



(48) 
(49) 
(50) 
(51) 
(52) 



The constant value at infinity and 1/r fall-off term in 
the aether are due to the unit timelike constraint. Thus, 
unlike ordinary fields, the aether will contribute to the 
energy expression directly. Inserting (|48)) - ([52l) into the 
von-Freud superpotential (|26p and aether superpotential 
yields the usual 'ADM mass" of GR 



Eqr — 



and the aether modification 



Ci4 



8ttG 

Combining, we find 



du n,d b- 2 



p _ , Ci4 

i'tot - ^(1 - — ). 



(53) 



(54) 



(55) 



This shows that the aether contribution effectively renor- 
malizes the ro/2G value we usually find for the total 
energy of an asymptotically flat spacetime in GR. This 
renormalization can also be understood as a rescaling 
of Newton's constant of the form Gn — G/{\ — C14/2), 
which agrees with the result of [22] . 

Equation (|55p implies that if C14 < 2 then the total en- 
ergy of the Einstein- Aether theory is positive if the ADM 
mass tqI^IG ispositive. However, the positive energy the- 
orem for GR 3 requires a stress-tensor that satisfies the 
dominant energy condition. The aether stress-tensor Q 
does not appear to generally satisfy this condition, so 
proof of total positive energy remains elusive. For some 
speculative thoughts on modifying the positive energy 
theorem, see Q. 



Despite these difficulties, there are special cases of the 
non-linear theory that are simple enough for calculations 
of the energy, yet still give important results. One sector 
of interest is the non- linear decoupled limit. As discussed 
above in Section |V] this limit allows one to essentially 
replace gab with the flat Minkowski metric r]ah in the 
aether parts of (HJ. One significant example is C2 = C3 = 
C4 = 0, ci 7^ theory. In this case the Lagrangian density 
for the aether is 



(56) 



This corresponds to a nonlinear sigma model on the unit 
hyperboloid, which has a stress tensor satisfying the dom- 
inant energy condition. A simple way to see this is to note 
that the derivatives of the individual scalar components 
and are contracted with 77^^, which is positive defi- 
nite on the unit hyperboloid. Returning to the linearized 
plane wave energy densities of Section |V] we see that in 
this special case of (p9)) - (|4T|) . if < ci < 1 energy is 
positive in both the linearized and decoupled non-linear 
regimes of the theory. 

Another important application of the decoupling limit 
relevant for our analysis of energy is the work of Clay- 
ton [2^. Clayton examined the Maxwell-like simplified 
theory where ci — —c^, C2 = C4 = in the decoupled 
version of non-linear Lagrangian ©, yield mg 



L= d 



!)}• 



, (57) 

where F^^ = — d^u^. The standard calculation of 
the Hamiltonian and the constraint equations then pro- 
duces the following on-shcU value for the Hamiltonian 



H= / d'xaP^ + P'd,uo + \{F,,f}. 



(58) 



Unlike the electromagnetic case, the second term cannot 
be turned into a total divergence since now V - P = — Aug 
on-shell. This implies that for some solutions the value 
of the Hamiltonian is negative. For example, as initial 
data choose Ui to be the gradient of a scalar field and 
Pi = —diUQ. Evaluating the Hamiltonian then yields 



E : 



-l/2(a,;Uo)2 



(59) 



which can be made arbitrarily negative by an appropriate 
choice of uo ■ 

Moreover, as Clayton points out, the negative energies 
are not restricted to this special case. In particular, al- 
lowing C2 7^ does not affect the P ■ d uq term in the 
Hamiltonian and even produces additional questionable 
terms. The indefinite nature of the decoupled Hamil- 
tonian contrasts with with the wave energy densities of 
Section |Vl which clearly can be made positive definite in 
the Maxwell-like case. The key point is that the wave 
results are in the linearized theory and associated with 
quadratic parts of the Hamiltonian, while the indefinite 
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terms appear at higher orders. For example, the hn- 

earized constraint equation v'^ = eliminates P ■ d uq 
from the Maxwell-like Hamiltonian and forces the uq in 
(|59p to be quadratic or higher in the perturbation. Thus, 
the indefinite pieces begin to appear at quartic order in 
the Hamiltonian. This indefiniteness at higher orders im- 
plies that the decoupled, linearized results of Lim and the 
"coupled" , linearized analysis of this paper generally do 
not detect possible energies of arbitrary sign in the fully 
non-linear decoupled Einstein- Aether theory. 



VII. DISCUSSION 

In this paper we have derived two energy-momentum 
pseudotensor expressions for the Einstein- Aether theory 
and used them to compute the energy densities of weak 
gravitational waves and the total energy of an asymptot- 
ically flat solution. The constraints of Section |V] show 
that a sector of this LV model satisfies the important 
theoretical condition of positive energy in the linearized 
case. However, a remaining open question is whether 
the energy remains positive when we consider the full 
non-linear theory. We have argued that in the decou- 
pled limit the ci ^ non-linear sigma model is immune 
to the sickness of energies of indefinite sign. However, 
other special cases of the coupling constants yield nega- 
tive energy solutions even when the linearized theory has 
positive energy. A complete answer to the question of 
positivity of energy in the non-linear theory is not yet in 
hand. 

The recipe for the Weinberg pseudotensor discussed 
in Section IHII and the Einstein superpotential and pseu- 
dotensor derived in Section IIVI also has applications in 
studying the emission of gravitational-aether radiation 
from astrophysical sources. In [2l| the analog of the 
quadrupole formula (which also involves monopole and 
dipole moments), is obtained using a pseudotensor ex- 
pression derived from the related Noether charge ap- 
proach. This expression is then used to track radiative 
energy losses and study constraints on the model. 
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ism of Julia and Silva [l8[ applied to Lagrangians that 
depend on the fields and their first and second deriva- 
tives. We can write a variation in the Lagrangian as 



5L = ^5^P 



dL 



dip d{d^i!) 



dL 

^^^^^^ + Wm^^^'^""^^ 



and then integrate by parts to isolate the equations of 
motion E and a symplectic current 9^ , 



5L ESiP + d^e'', 



(A.2) 



where 



dL 



dL 



d{df,d„tp) 



(A.3) 



If the action associated with L is invariant under a con- 
tinuous transformation of the fields. 5L = d^S^. Thus, 
we have the equation 



5^(5'^ -61") = E5i!. 



(A.4) 



This identifies the on-shell {E = 0) conserved Noether 
current, 



dL 



d{d^dui}}) 



d{d,rp) 
d,{SiP)~S^. 



'd{d^duip)' 



(A.5) 



We now want to consider a gauge transformation of the 
fields that involves derivatives of the generator S,^{x). 
Here we will focus on the special case restricting atten- 
tion to only the first derivative. Following the analysis 
and notation of plsf we parameterize the gauge transfor- 
mation as 

6^ = ^^Aa + id,i^)A''A (A.6) 



where A is an internal or spacetime index and A is 
a transformation matrix. The quantity can be ex- 
pressed similarly as 



Appendix 
Background 



5" = e-^E^; + (9.e'^)E7 + idrd,^^)^''/''l (A.7) 



Inserting these forms into (jA.Sp and combining terms, we 
find that on-shcU 



In this appendix we will derive the Einstein psuedoten- 
sor and superpotential using the Noether current formal- 



d, {^^j^ + d.i^u^' + drd.ev^^''^^) = 0, (A., 
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where 



dL 



tI.(Tu) 



dL 



A, 



(A.9) 



dL 



dL 



d{dpdrip) 



dL 



d{d^,dri)) 



fJ.{Tu) 



-Aa 
(A.IO) 
(A.ll) 



Since and its derivatives should be arbitrary and in- 
dependent, this single equation decomposes into 4 equa- 
tions 



a scalar density C and a total divergence d^j^W^, 

C{ip, dip, aV) = di!, 5V) + ^^,[W^'{■^p, dip)]. 

(A.19) 

If is a vector density then the total divergence is 
a scalar density, but we allow for a non-covariant to- 
tal divergence. For a variation that is an infinitesi- 
mal diffeomorphism generated by a vector field , we 
have 6C — dp{^^C) since £ is a scalar density and 
5{dpW^') = dp[5W^'). Therefore the surface term 5^ 
in IA.4I has the form 



Now consider the Einstein- Hilbert action 



5'eh — 



-gR 



(A.20) 



(A.21) 





« 


(A.12) 


J% + d^Uj 


« 


(A.13) 


u^^"^ + drVl^^"^ 


= 


(A.14) 


^ A 


= 0. 


(A.15) 



The first two equations hold on-shell, while the last two 
are identities (since there are no second or third deriva- 
tives of on the right hand side of (|A.3[) V 

The gauge symmetry implies that is conserved and 
equal to the divergence of the superpotential . Since 
= ^"^(2;), the Noether current J*^ (jA.Sp will now be 
parameter dependent. Let us consider a one parameter 
subgroup of the local gauge or diffeomorphism symmetry 
where (/^ has the decomposition. 



(A.16) 



and is fixed. Inserting this form into (jA.13|l produces 



with J^^ = ^t^^- The conserved charge is 



Q 



J? d^x ■ 



?0 ^ . 



(A.18) 



Q depends on the choice of and can be expressed in 
terms of the gauge fields up using (|A.10[) . If is an 
asymptotic translation in an asymptotically flat space- 
time, then the conserved charge will be a total energy or 
momentum. Using the variational principle {6S — ^ 
equations of motion), we will show in the next section 
that the choice of S,q is subject to certain boundary con- 
ditions at infinity. We now apply this type of Lagrangian 
analysis to the Einstein- Aether theory. 



of pure GR. The Ricci scalar R has a dependence on 
second derivatives of the metric. In light of this, Einstein 
exploited a property of the Hilbert action that allows it 
to be separated into a bulk and a surface term 



R d X — I \/—gL 



^bulk 



d„y^'d'^x. 



(A.22) 



This decomposition of the Ricci scalar takes the following 
form 



ibulk 



~ 9 aS^ r)p 



^ »7<5^ a/3/ 



(A.23) 



~9{Kp9"^ -^ipgn- (A.24) 



where T is the Levi-Civita connection. One can elimi- 
nate the total divergence by adding its negative to the 
Einstein-Hilbcrt action 



y V^Lbuik = y i? - a^y'^d^^. (A.25) 



The elimination does not affect the equations of motion 
and is consistent with the general action (jA.19p with C = 
y/—g R and — —V^. The result of this is a loss of 
diffeomorphism invariance since the remaining Lbuik in 
the 'T^" action is not a scalar. We have allowed for this 
possibility with the non-covariant d^W^ term in (|A.19p . 

With C{ip,dip,d'^ip) the bulk part of the Einstein- 
Hilbcrt action plus the the aether terms in ([T]), we arrive 
at the Einstein- Aether form of (|A.4p on shell 



da 



dC 



d{df,gaf3 



■Sgal3 



dL 



d{dpu- 



-6u" - S^' 



0. (A.26) 



Application 

Assume that the Lagrangian density C{il!,dip,d'^ip) is 
invariant under diffeomorphisms and is a combination of 



The second derivative terms in (jA.5|) vanish in this case. 
Under a diffeomorphism generated by a vector field 
the variation of the metric and the aether is simply the 
Lie derivative, dgaf3 = i^d^g^p + daS,''g„p + dp^^„ and 
,5m" = C'd^u"' - u-^dj^". It follows from (fOOl) that 
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S'' = ^^Uuik - V^idsd^^g^" - dsd^i'^g'^^). Inserting 
these forms into (|A.26p produces 



f 



dL 



d„gaf3 



dL 



dL 



d{d^ga0) 



dL 



(A.27) 



(A.28) 



as coefficients of d£,^ and d^£,^ respectively. t,j^ and 
C/^^^, and K ^(t^) are the analogs of J^, t/j^'', T/j'^'^^) in 
(|A.9|l - (|A.lip . The resulting equations due to the arbi- 
trariness and independence of the derivatives of ^'^ are 



showing that a Noether charge is again obtained as a 
surface term. Einstein effectively chose the vector to 
be a constant in (jA.34p . reducing the pseudotensor to a 
form consistent with the flat spacetime canonical stress 
tensor (fTg]) . However, this choice is not inconsequential. 
The variational principle for the action requires the 
vanishing of the surface term in the asymptotic region, 



dC 



didf,gaf3 



■5gai3 



dC 



did^W-) 



(A.35) 



Therefore we have Dirichlet boundary conditions Sga/j = 
and = on the metric and the aether at infinity. 
Inserting the Lie derivatives for the variations above, we 
see that as r — + oo 



d t 
f 







(A.30) 
(A.31) 
(A.32) 
(A.33) 



Following (jA.16|l . we can keep the vector fixed (and 
determine it later for each conserved charge) by choosing 
^'^ = e(x)^g. The main result, as before, is 



(A.34) 



^ 
0. 



(A.36) 
(A.37) 



Since f has been chosen to be constant everywhere and 
is asymptotically constant, the connection coefficients 
must vanish as one approaches spatial infinity. Thus, 
one must compute the pseudotensor and superpotential 
in a coordinate system where the connection vanishes 
asymptotically as 0(l/r) or faster. 
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